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We demonstrate that the strong anomalies in the high frequency LO-phonon spectrum in cuprate 
superconductors can in principle be explained by the enhanced electronic polarizability associated 
with the self-organized one dimensionality of metallic stripes. Contrary to the current interpretation 
in terms of transversal stripe fluctuations, the anomaly should occur at momenta parallel to the 
stripes. The doping dependence of the anomaly is naturally explained, and we predict that the 
phonon line-width and the spread of the anomaly in the transverse momentum decrease with in- 
creasing temperature while high resolution measurements should reveal a characteristic substructure 
to the anomaly. 



The phonon spectrum of the high-T c superconducting 
cuprates is characterized by a peculiar anomaly: halfway 
the Brillioun zone the 'half-breathing' Cu-0 vibration 
mode seems to suddenly dip down to a much lower fre- 
quency [3, 0, 0] • The line- width reaches its maximum at 
a wave-vector somewhat shifted from the frequency-dip 
position, while it narrows at higher temperatures [4|]. In 
addition, the anomaly has a narrow intrinsic peak width 
as function of momentum transversal to the mode prop- 
agation direction [4] . This is hard to explain in a conven- 
tional fermiology framework, and interpretations invok- 
ing a coupling between the phonon and purely electronic 
collective modes of the stripes acquired credibility [5[ by 
the recent demonstration that the anomaly is particu- 
larly pronounced in La 2 -i/sBai/ 8 Cu04, [3U4j], a system 
with a well developed static stripe phase [6(. Initially it 
was believed that the relevant stripes modes should corre- 
spond with the transversal 'meandering' fluctuations 0, 
a possibility further quantified by Machida et a/.0 by 
computing the Gaussian fluctuations around the Hartree- 
Fock stripe ground state. This interpretation has how- 
ever run into a serious objection: the phonon anomaly 
is also observed in the YBCO superconductor, where the 
orientation of the stripes relative to the orthorombic lat- 
tice can be directly deduced from the anisotropy of the 
spin-fluctuations in the untwinned crystal. It turns out 
that the phonon anomaly occurs for phonon- wavevectors 
parallel to the stripes, at a right angle as compared to 
the expectations for transversal stripe modes 0. 

A well known problem with the Hartree-Fock stripes is 
that they are 'insulating like' jlOj, while cuprate stripes 
are quite metallic. This supports the idea that they form 
an electronic liquid crystal of the smectic kind as intro- 
duced by Kivelson and coworkers pjj . Here the transver- 
sal modes are frozen out by commensuration effects and 
instead the low energy physics is governed by on-stripe 
compressional fluctuations. As one of us already noticed 
some time ago [12], the electronic polarizability associated 
with the on-stripe Luttinger-liquid like physics should be 



strongly enhanced at the on-stripe 2k p wave- vectors and 
this feature might well govern the phonon anomaly (see 
Fig. [Tj). Embedding this ID electron dynamics in a 2D 
optical phonon background is still a non-trivial exercise 
and here we present a first attempt to compute these mat- 
ters quantitatively. We employ two simplifying assump- 
tions: (i) we use the free fermion charge susceptibility 
(Lindhardt function) instead of the fully interacting 'Lut- 
tinger liquid' form, since the phonon anomaly appears to 
be in first instance sensitive only to the gross features of 
the ID electron dynamics, (ii) More critical, we assume 
that at the energy of the anomaly the inter-stripe interac- 
tion effects have already diminished and this might well 
be an oversimplification. The outcomes are as follows. 
The stripe alignment problem is solved by construction, 
but we also arrive at some less obvious predictions: (a) 
the anomaly is now caused by the phonon crossing the 
ubiquitous continuum of ID charge excitations centered 
at the intrastripe 2k p (Fig. [T]). At the phonon frequency 
uj this continuum has a momentum width Aq w kpuj/v c 
where v c is the electronic charge velocity and this causes 
a 'double dip' structure in the phonon spectral function 
(Fig. 2) , not unlike to what is predicted by Machida et 
at. [8[ (inset Fig. [2]). However, the big difference with 
the propagating transversal modes of Machida is that in 
the 'smectic scenario' the phonon is completely (Landau) 
damped in the momentum region in between the dips: 
high resolution neutron scattering measurements should 
be able to resolve this, (b) Since the characteristic mo- 
mentum where the anomaly occurs is now determined by 
the intrastripe electron density and not by the interstripe 
distance the position of the anomaly should have a dop- 
ing dependence that is radically different from what is 
expected for transversal modes (Fig. [3]). (c) We predict 
that the anomalous phonon linewidth Auo decreases with 
increasing temperature (FigHj).(d) Finally, although we 
assume that the stripe Luttinger liquids are strictly inde- 
pendent, the form factors of the electron phonon interac- 
tion [13] localize the phonon anomaly in 2D momentum 
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space also in the direction perpendicular to the stripes 
(Fig. [4]). Counterintuitively, we find that it contracts 
further when temperature is raised (inset Fig. a be- 
havior that is also observed in experiment [4| , while it is 
difficult to rationalize invoking transversal stripe modes. 

We embed the array of parallel metallic stripes in the 
2D phonon universe by considering a simplified propaga- 
tor: it describes non-interacting electrons in Bloch states 
on a periodic array of parallel lines in two dimensions, 
moving freely along these lines, but with vanishing inter- 
line overlap of the 'Wannier' functions: 



k M'=o 



N ~ 1 ik(f-r')+iQ(ly-l'y') 



iuo — e(k) 



(i) 



where r, r' span the the ID Cix-centered stripes(lines) 
in real space; uo is the Matsubara frequency. We con- 
sider an orthorombic Cu-0 plane, with a(b) being the 
unit cell spacing along the x(y)-axis. The stripe metal- 
lic direction is along the x-axis and the stripe Umklapp 
momentum Q = 2n/bN is along the ?/-axis, where N 
(= 4 at higher doping) is the number of unit cells in one 
inter-stripe (charge-density) period. The momentum k 
spans TV/e sites in the reduced orthorombic Brillouin zone 
< k y < 2ir/bN; < k x < 2ir/a. The electron dis- 
persion e(k) ~ e(k x ), ignoring inter-stripe tunneling (see 
below) . 

A hole on site r inside the stripe communicates 
with neighboring oxygens bond-stretching displacements 
u r ±i = Ui(r =b i/2), i = a, b and we take the effective 
coupling between those and the Zh ang -Rice singlets as 
introduced by Horsch an Khaliullin [131 ]. 



H e - ph = g ( U x ~ u -x + U y~ U -y) C W ( 2 ) 



with go w 2eV/A, using the standard estimates for the 
charge transfer energy and hoppings[3]. This is actually 
our main step: from this Hamiltonian and 2D 'striped' 
electron propagator Eq. (pQ) it is straightforward to cal- 
culate the self-energy part of the dynamical matrix asso- 
ciated with the Cu-0 plane: 



A^(q,u;)=uooU(q x ,u) 



(3) 



where uoo is bare phonon frequency, s qx = sin^a/2, 
s qy = smq y a/2, and a,/3 = 1,2,3 enumerate ions in 
the in-plane Cu-0 unit cell: this particular form of the 
electron-phonon form-factors follows immediately from 
the tight binding Hamiltonian Eq. ([2j), implying actu- 
ally a substantial dependence on the momentum q y per- 
pendicular to the stripes. H(q x ,u;) corresponds with the 
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FIG. 1: Real- (black line) and imaginary (red line) part of the 
on-stripe electronic phonon self-energy U(q x ,cu) at T = 10K 
at a fixed (phonon) frequency u> — 68meV as function of the 
momentum component q x parallel to the stripes, taking for 
the dimensionless electron-phonon coupling a value represen- 
tative for the cuprates: £ = 0.5. In the inset the results are 
shown at a higher T = 300K. 



polarization propagator of the (in principle, interacting) 
on-stripe Luttinger liquid, depending on the momentum 
component q x along the stripes, while its fermion lines 
are given by the propagator Eq.([T]). 

This self-energy has to be added to the bare (undoped) 
afr-plane ionic 6x6 dynamic matrix A^ ,/3 (<f, u) of the or- 
thorombic planes, constructed to be in close agreement 
with the experimental data Q for the in-plane bond- 
stretching LO phonon modes in the undoped cuprates. 
The phonon spectra w a (q) associated with polarizations 
are obtained by the diagonalization of the total dy- 
namic matrix [15[. From these solutions the phonon 
spectral functions given by the imaginary part of the 
phonon propagator D(q,cu) are obtained. Notice that 
the A^'s are in general complex quantities, with the ef- 
fect that oo a (q) acquires an imaginary part representing 
the phonon damping. 

We use for H(q x ,w) the well known ID Lindhardt 
function lldl, 
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here A± = (2p ± q) 2 q 2 — oj 2 , all momenta and energies 
are measured in units of the Fermi-gas parameters kF and 
6f = k F /2m ~ leV; r = ksT/cF is a dimensionless tem- 
perature. The dimensionless electron-phonon coupling 
constant £ = g^/KcF ~ 0.5 is representative for cuprates 
with K w 25eV/ 'A 2 the lattice force constant [Uj]. 

The effect of this ID polarizability on the phonons 
follows from the behavior of H(q x ,u) at the phonon- 
frequency uo = uoq as function of q x (Fig. 1). The contin- 
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FIG. 2: False color plot of the LO phonon spectral function 
vs momentum in the stripe direction and energy at low tem- 
perature (10K) for ID 6f ~ leV and the same parameters as 
in Fig. [1] Phonons couple merely to the ID Fermi-gas like 
charge excitations of the metallic electron system confined in 
the stripes. Strong phonon damping in a momentum region 
'inside' the anomaly is caused by the decay in the continuum 
of quasi ID electron- hole excitations. Left inset: different be- 
havior of the 'standard' mode coupling of the phonon with 
a propagating stripe collective mode (e.g., ref.0). Right in- 
set: the renormalized phonon dispersions determined by the 
crossing of the phonon frequency ujo and the real part of the 
ID polarization propagator of Fig. [T] see text. 



imm of charge excitations in a ID fermi-gas is the well 
known fan in the momentum-frequency plane, centered 
at 2kp at uo = and bounded by 2kp ± uo/vp at fi- 
nite frequency. For non-interacting electrons the spectral 
function (/mil) is just the box of Fig. (pQ), while in the 
presence of interactions the spectral weight will pile up 
at the edges. Since II is proportional to the phonon self- 
energy, the phonon spectral function, Fig. ([2|), indicates 
that the phonon dispersion is pushed downwards when it 
approaches the edges of the quasi ID electron-hole con- 
tinuum from either side, to broaden strongly when it en- 
ters the continuum. This is markedly different from the 
result based on a mode coupling between the phonon and 
a propagating mode as for instance discussed by Machida 
et al. [8[ as shown in the inset of Fig. ([2j): in this 
case there is no phonon damping and the intensity is 
just distributed over the propagating modes subjected to 
an avoided level crossing. 

The results shown in Fig. ([2]) suggest that it is unprob- 
lematic to give a quantitative account of the observed 
electron-phonon effects in terms of the metallic stripe 
electrons, provided that flatness of their band dispersion 
e ±(ky) perpendicular to stripes is bounded from above 
by ~ uoo. The numerical LDA+Z7 calculations [17] sup- 
port this assumption, suggesting e± ~ lbmeV « 0.2uo. 
Altogether, this success is not all that remarkable since 
the 'propagating mode' and the ID metallic stripe sce- 



nario's give in this regard the same answer: by concen- 
trating the charge excitations in a small kinematic region 
(as compared to conventional fermiology) one is dealing 
in essence with an avoided mode crossing story and be- 
cause this is about resonance small numbers can give big 
effects. 

A distiction between the effects of the transversal 
stripe modes and the internal ID-like fermionic excita- 
tions on the phonon anomaly should be revealed by the 
different doping dependences of the locus of the phonon 
anomaly in momentum space. The transversal fluctu- 
ations emerge at the stripe ordering wave vectors and 
these should follow the famous Yamada plot[18[, corre- 
lating the stripe ordering wavevectors 5 with doping x 
(Fig. [3]), such that at low dopings the anomaly should 
live at a wavevector q a ~ \ jx. On the other hand, deal- 
ing with the quasi ID modes the locus of the anomaly is 
determined by the on-stripe electron density and accord- 
ing to the Yamada plot this stays constant ('half- filled') 
at q a = 2kp = 7r/2a up to x c = 1/8, while at higher 
dopings it should follow q a = 2kp oc x — x c because the 
on stripe hole density is increasing (the blue line in Fig. 
[3]). Experimentally the locus of the anomaly in /c-space 
is conspicuously doping independent [3[ , actually arguing 
strongly against the transversal mode. It would be inter- 
esting to find out if the 'center of mass' of the anomaly 
does shift at higher dopings. 

Our theory yields a rational for the observed gross tem- 
perature dependences of the anomaly [4]. The rather 
counterintuitive narrowing of the frequency width with 
increasing temperature follows naturally from the ID po- 
larization propagator decrease, particularly: /mil oc uo/T 
(Fig. [I]). The effect of this change on the phonon-spectral 
function is shown in Fig. (|4]): a substantial narrowing oc- 
curs at higher temperature. Moreover, right at the 'dips' 
a substantial phonon hardening occurs since the phonon 
positions at these momenta are most sensitive to the de- 
tails of the real part of the self energy. 
The phonon anomaly behavior in the 'transversal' q y - 
direction (inset, Fig. [4j) follows from our analysis 
of the expression for the phonon spectral function 
ImD(q x , q y \uj) showing a substantial q y dependence close 
to the anomaly due to the form-factors in Eq. (|3|): 



ImD(q x ,q y ;uo) 
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and uo, 



q x ,q y 



U0 



0.5A qx:qy ReH qx is the renormalized optical phonon mode 
frequency. The width of the Lorentzian with respect to 
q y at uo = oj qx =2k F ,q y =o is: Sq y « 2/a y / \ImU/ ReU\ ~ 
0.1 x 27r/a, assuming a fiat bare mode dispersion. Given 
the ratio under the square root, the width 5q y decreases 
when the temperature is raised. 

In summary, we have analyzed a minimal- and likely 
too oversimplifed model dealing with the 'Luttinger 
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FIG. 3: The qualitatively different doping dependences ex- 
pected for the characteristic wave vector of the anomaly when 
transversal stripe model (1) or the intra- stripe quasi ID elec- 
tron excitations (2) are responsible. In case (1) the anomaly 
should follow the stripe-ordering incommensurate wave vec- 
tor 5 dependence on doping x: the famous 'Yamada plot' [3 
(black dots) ; (blue line) is for the softened phonon wave- vector 
q = 2Uf in case (2). 
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FIG. 4: Calculated phonon spectral function ImD(q x ,q y — 
0]uj) dependence on frequency cj, using the same parame- 
ters as in Fig|2j at different temperatures. A temperature 
narrowing is predicted, accompanied by a phonon hardening 
at the 'dips' seen in Fig(2] Inset: calculated ImD(q x ,q y ;cu) 
momentum dependence in the 'transversal' q y direction per- 
pendicular to the stripes, at fixed q x = 2k f and co = 68 meV. 
Counter intuitively, the anomaly even localizes further when 
temperature is raised, in accord with experiment 0]. 



liquid'-like excitations coming from the electrons con- 
fined in stripes interacting with optical lattice phonons. 
The gross features of the phonon anomaly as measured 
experimentally are consistent with the workings of a 
quasi ID array of metallic intrastripe Luttinger liquids. 



More importantly, although it will take a significant in- 
vestment, it does make sense to measure the phonon 
anomaly with a much higher resolution than has accom- 
plished up to now. Such measurements should reveal the 
internal structure of the anomaly and as we showed in 
this paper, from this information much can be learned 
regarding the still quite mysterious organization of the 
electron systems in high T c superconductors. 
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